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Relativistic chaos is coordinate invariant 
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The noninvariance of Lyapunov exponents in general relativity has led to the conclusion that chaos depends 
on the choice of the space-ftme coordinates. Strikingly, we uncover the transformation laws of Lyapunov expo- 
nents under general space-time transformations and we find that chaos, as characterized by positive Lyapunov 
exponents, is coordinate invariant. As a result, the previous conclusion regarding the noninvariance of chaos 
in cosmology, a major claim about chaos in general relativity, necessarily involves the violation of hypotheses 
required for a proper definition of the Lyapunov exponents. 
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PACS numbers: 98.80.Jk, 04.20.-q, 05.45.-a 

Chaotic properties of dynamical systems with a 
reparametrizable time coordinate are important in phys- 
ical theories without an absolute time, such as general 
relativity. The study of chaos in general relativity has 
followed two main lines. One considers the geodesic motion 
of test particles in a given gravitational field [ 1 ] . The other 
investigates the time evolution of the gravitational field itself 
001, which is relevant in cosmology. While the former case 
has been studied with standard methods of the dynamical 
systems theory, an adequate characterization of chaos in the 
latter is currently an open problem. The difficulty comes 
from the dependence of the shear between nearby trajectories 
on their time parametrization. Accordingly, dynamical 
properties, such as mixing and initial-condition sensitivity, 
may depend on the time parametrization. 

In classical physics, the study of dynamical systems con- 
cerns differential equations of the form 



dx/dt = F(x), 



(1) 



where t is a uniquely defined parameter that usually represents 
the time. Although a general definition of chaos is missing, it 
is widely accepted that it regards the dynamics of bounded 
orbits and that a chaotic system must present sensitive depen- 
dence on initial conditions @]. Chaos can then be quantified 
in terms of Lyapunov exponents insofar as the following 
conditions are satisfied: (z) the system is autonomous; (ii) the 
relevant part of the phase space is bounded; (Hi) the invari- 
ant measure is normalizable; (z'v) the domain of the time pa- 
rameter is infinite. Such a characterization is convenient be- 
cause it is invariant under space diffeomorphisms of the form 
y = i/>(x). As a result, chaos is a property of the physical sys- 
tem and does not depend on the coordinates used to describe 
the system. 

In general relativity, the nonexistence of an absolute time 
parameter forces us to consider Eq. {0 under space-time dif- 
feomorphisms: y = tj}(x,t), dr = X(x,t)dt. A conceptual 
problem then arises because of the dependence of classical in- 
dicators of chaos, like Lyapunov exponents and entropies, on 
the choice of the time parameter. This problem has attracted a 
great deal of attention since it was first identified in the mix- 



master cosmological model ((J, whose largest Lyapunov ex- 
ponent was shown to be positive or zero for different choices 
of the coordinates |7]. In particular, numerous methods based 
on invariant curvature |8], symbolic dynamics |9], Painleve 
analysis ifioll . and fractals 

OH 

, have been proposed to- 
ward an invariant characterization of chaos in cosmology. The 
problem, however, goes beyond relativistic cosmology since 
it has been argued that the same kind of noninvariance can 
be observed in a system as simple as the harmonic oscilla- 
tor if time reparametrizations are allowed lll3ll . Moreover, it 
has been exhibited examples of systems whose nonmixing dy- 
namics can be converted into a mixing one through a time 
reparametrization 11411 . These results have led to the tacit as- 
sumption that chaos itself depends on the space-time coordi- 
nates. In general relativity, this noninvariance would imply 
that chaos is a property of the coordinate system rather than a 
property of the physical system (see Ref. |3|. and references 
therein). 

In this Letter we investigate the transformation laws of 
Lyapunov exponents under space-time reparametrizations that 
preserve conditions (z'-z'v). To be specific we consider an Eu- 
clidean phase space (not to be confused with the pseudo- 
Riemannian spacetime), where the relevant invariant measure 
is the natural measure. Our principal result is that Lyapunov 
exponents transform according to 



K/(\)t 



(i = l,...,N), 



(2) 



where < (X)t < 00 is the time average of A = dr/dt over 
typical trajectories and N is the phase-space dimension fl5ll . 
The values of the Lyapunov exponents are, of course, nonin- 
variant because a simple reparametrization such as (x, t) — > 
(x, at) transforms the exponents h l into K l jot. However, the 
signs of the Lyapunov exponents are invariant. In particular, if 
h > is the largest Lyapunov exponent, the reparametrization 
only changes the characteristic time-scale T = 1/h for the 
manifestation of the chaotic behavior. The striking implica- 
tion of our findings is that chaos, as characterized by positive 
Lyapunov exponents, is coordinate independent. As we show, 
the vanishing of Lyapunov exponents in the mixmaster cos- 
mology as well as in other examples previously considered 



2 



in the literature is due to the violation of at least one of the 
conditions (z'-z'v) above, which are required for the Lyapunov 
exponents to be meaningful indicators of initial-condition sen- 
sitivity and hence chaos. For instance, it has been frequently 
claimed that the exponential divergence of initially close tra- 
jectories that separate as Sx(t) = 6xoexp(ht) can be re- 
moved with a logarithmic reparametrization of the time 1 16]. 
The suggested transformation is defined as t = In t, so that 
5x(t) = 5x.oT h . This reparametrization is then interpreted 
as converting a positive Lyapunov exponent into a zero Lya- 
punov exponent, and the inverse of this transformation has 
been used to support the claim that integrable systems can 
have positive Lyapunov exponents 1 1 3 ] . The problem with 
this argument is that, if the original system is autonomous, the 
reparametrized system is necessarily nonautonomous. Alter- 
natively, if we increase the dimension of the phase space in 
order to eliminate the explicit dependence on time, the orbits 
become unbounded. In any case, Lyapunov exponents are not 
valid indicators of chaos. We show that a similar problem, al- 
though more subtle, is present in the mixmaster cosmological 
model. 

The Lyapunov exponents of an invariant set of the phase 
space of system Q are defined as 



t-oo t It/So I 



(i= 1,...,A0, (3) 



where dr}\(t)/dt = DF(x(f)) • r]\(t), 



(4) 



x(0) is a typical initial condition, and T]l — r\\{<S] are tangent 
vectors at x(0). We assume that F is a continuously differen- 
tiable function of x and system Q has N independent Lya- 
punov exponents. Behind definition (13141 are the hypotheses 
(i-iv), namely, that with respect to (x, t), the function F does 
not depend explicitly on t, the dynamics is well defined for 
t in the interval [0, oo), and the invariant set is bounded and 
has finite natural measure 1 17]. Otherwise, positive Lyapunov 
exponent is not a well defined criterion for chaos. To ensure 
that system Q remains autonomous after the coordinate trans- 
formation (x, t) — > (y, r), we consider that when functions 
A or ip depend explicitly on t, the coordinates x and y are 
redefined to incorporate t as an additional dimension in the 
phase space fliil . As a result, the coordinate transformation 
is always reduced to a time-independent transformation of the 
form: 



y = t/>(x), dr = X(x)dt, 



(5) 



where A is a strictly positive, continuously differentiable func- 
tion, and ij) is a diffeomorphism. This is the general class of 
transformations for which integrability is coordinate invari- 
ant III in the sense that, if I2, . . .} are independent in- 
tegrals of motion with respect to the coordinates (x, t), then 
{Ii o ip~ 1 , 1 2 o ip~ 1 , . . .} are independent integrals of motion 
with respect to (y, r). 

Transformation (|5j is composed of a time reparametriza- 
tion followed by a space diffeomorphism. It is well known 



that the Lyapunov exponents are invariant under space dif- 
feomorphisms |4]. Without loss of generality, we consider 
only transformations of the time parameter: (x, t) — > (x, r), 
where dr — X(x)dt. All the orbits of the phase space are 
invariant under this kind of transformation as the velocity 
field of the reparametrized flow is parallel to the original one: 
dx/dr — A _1 (x)F(x). The Lyapunov exponents, however, 
may be different because Eqs. J3l4i become 



K = lim Iin^l 
r^oor \r) l r0 \ 



(6) 



and dr,Ur)/dr = D[A- x F](x(r)) • tj* (r), (7) 

respectively. Incidentally, mixing, which is a property most 
often observed in chaotic systems, is not invariant under trans- 
formation since it has been shown that an adequate time 
reparametrization of a nonchaotic irrational flow on a 3 -torus 
is mixing IU4I1 . 

For the same initial conditions, the defining relations of h\ 
and h\ present two different factors. The first, associated with 
the time average of A, comes from the difference between the 
two time parametrizations along the same orbit and is factored 
out when Eq. is written as 
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\Vr(t)\ 

K \ 



(8) 



where = ^OK*)), r(t) = f* X(x(t))dt, and (A) t = 

liirit^oo j J Q X(x(t))dt. Condition < (X) t < 00 is a basic 
requirement for the natural measure to be well defined 11711 . 
The second factor, due to the gradient of A, is associated with 
the difference between the surfaces of simultaneous time for 
each parametrization and is separated out when rf r in Eq. @ 
is parametrized in terms of t rather than r: 



dt 



= DF (x(t)) ■rf i T (t)—[F-V i In A] (x(t)) • rf T (t) . (9) 



The last term in this equation implies that the time evolution 
of vector rf T (t) in Eq. (|8|i is in general different from that of 
vector r\\ (i) in Eq. (|3j. But the relevant question is: how large 
is this difference? 

Here we show that the difference is sub-exponential, in the 
sense that rf t (i) and rf T (i) grow or shrink with the same ex- 
ponential rate. This implies 



limimMWi 

t->oo t \r)l Q \ 



(10) 



which in turn implies our main result (|2}. 

First we analyze periodic orbits, which form the most fun- 
damental building blocks of chaotic sets 12011 . On a peri- 
odic orbit x*, it is convenient to adopt the explicit notation 
rjl(x*(t),t) and rf T (x*(t),t) for the solutions of the varia- 
tional Eqs. (0} and 0, respectively. If x* is a fixed point, 
we trivially have rf T (x* (t) ,t) — rf t (x* (t) , t) because the 
last term in Eq. (|9} is zero. Now consider that x* is a peri- 
odic orbit with least period T > with respect to t, so that 
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x*(t + T) = x*(t). Let h\(-x*) denote the local Lyapunov 
exponents for Eq. Q on x*, and crj(x*) = exp(/ij(x*)) de- 
note the corresponding local Lyapunov numbers, where i = 
1, . . . , N. One of the Lyapunov numbers, say , is 1 because 
the bounded function dx* jdt is a solution of Eq. @, render- 
ing zero to the corresponding Lyapunov exponent. The same 
is true for any parametrization. To study the other Lyapunov 
numbers, let ix be the hyperplane orthogonal to F(x*(0)) at 
x*(0), and M : U C it i— > tt be the first return map on this 
hyperplane, defined in a neighborhood U of x*(0). This map 
does not depend on the time parametrization of the contin- 
uous flow, being exactly the same for both the original and 
reparametrized flow. It follows from standard results in Flo- 
quet theory |21] that the local Lyapunov numbers a % M (x*(0)) 
of this map, defined as the magnitude of the eigenvalues of 
the Jacobian matrix of M at x* (0), are the power T of the first 
N — 1 local Lyapunov numbers of the flow, i.e., a l M (x* (0)) = 
crJ(x*) T for i — l,...,N — 1. Geometrically, the local 
Lyapunov numbers of M can be defined as a l M (x*(Q)) = 
exp(h' l M (x*)T), where h % M {x*) is defined through Eq. (J3 
with r)l(x*(t),t) replaced with its orthogonal component to 
F(x*(t)). From the identity /i J M (x*) = hj(x*) then follows 

h l M (x*) = limt-.oo | In The same is true for 

r/*(x*(t), t), so that if h\ denotes the limit in Eq. (|8jl then 
?4(x*) = ~h\{x*) and ^ M (x*) = lim^oo \ In ■ 
From these, it follows that Eq. (1101 holds on periodic orbits. 

But the same must be valid in general because the Lyapunov 
exponents of typical orbits are weighted averages over the lo- 
cal Lyapunov exponents of all periodic orbits in the respec- 
tive ergodic component. The weight is the fraction of time 
spent by a typical orbit near the corresponding periodic orbit 
and is uniquely determined by the largest local Lyapunov ex- 
ponent (assumed to be positive). The smaller this Lyapunov 
exponent the longer it takes for the orbit to move away from 
that neighborhood. In other words, the largest local Lyapunov 
exponents define a natural measure, over which all the Lya- 
punov exponents are computed. Since h\{x*) — h\{x*), both 
the measure and the local exponents are the same for h\ and 
h\. Therefore h\ = h\ on typical orbits, and this is equivalent 
to Eq. E5- 

We now show that similar arguments can be extended di- 
rectly to typical orbits. Eq. @ can be written as a map 
M(r,t(t)) = nl(t + St), where r,\{t + St) = r,\(t) + 
<5tDF(x(i)) • r)\(t). For the natural measure to be normal- 
izable, F(x(t)) and V|F(x(f))| must not grow exponentially 
and F(x(i)) must not shrink exponentially along typical or- 
bits. Since dx/dt is a solution of Eq. (|4}, the component 
P(t]1) of r)\ parallel to the flow remains parallel in each itera- 
tion, and the Lyapunov exponent in this direction is zero. The 
orthogonal component Q{r}\) of rf t is mapped into two parts, 
one parallel and the other orthogonal to the flow. The paral- 
lel part is PMQ{rf t {t)) = St[F ■ V f In |F|](x(i)) • Q{r}\(t)), 
where neither F(x(t)) nor V In |F(x(f )) grows exponentially 
with t. If h\ < (hi > 0), each of the projections P(rf t ) and 
Q(f}\) must shrink at least (grow at most) as exp(ft,|i). But 



Q{r)\) cannot shrink faster (grow slower) than P(rjjf} because 
otherwise the contribution of the orthogonal part to Eq. Q 
would be negligible, resulting in h\ = 0, which violates the 
hypothesis that h\ ^ 0. Therefore, Q(r] l t ) ~ exp(/i|t) for 
both h\ < and h\ > 0. A similar relation is valid for if T (t) 
because \~ 1 dx./dt is a solution of Eq. (|9j. We then compare 
the solutions of Eq. (|9jl with those of Eq. for identical ini- 
tial conditions. The term involving matrix — [F- V 1 " In A] (x(t)) 
is parallel to the flow and does not affect the orthogonal part 
in Eq. l|9}. The term that contributes to the orthogonal com- 
ponent, DF(x(i)), is the same as that in Eq. 0. Therefore, 
Q(r)\(t)) = Q(rii(t)), which again leads to Eq. 

An important implication of our findings is that positive 
Lyapunov exponents are necessarily mapped into positive 
Lyapunov exponents under time reparametrizations. This im- 
plies that the previous examples of noninvariant chaos in cos- 
mology are based on the violation of hypotheses required 
for an interpretable computation of the Lyapunov exponents. 
Consider, for example, the mixmaster cosmological model 
la], which is believed to describe generic cosmological sin- 
gularities, and whose relevant hypotheses can be discussed 
explicitly. In the asymptotic limit (close to the big bang), 
the essential features of the continuous dynamics are repre- 
sented in the Farey map, F(u) = u — 1 if U > 1 and 
F(u) = u^ 1 — 1 if u < 1, whose Lyapunov exponent is 
zero. This result is claimed to be in conflict with the corre- 
sponding result for the first return map on [0, 1] (Gauss map), 
G(v) = 1/v — [1/v], where [l/v] is the integer part of 1/v, 
whose Lyapunov exponent is positive H: /i = 7r 2 /61n2. The 
problem here is that, different from the first return map, the 
orbits of map F are typically unbounded. Map F can be com- 
pactified for w = (u + by defining H : [0, 1] -> [0, 1], 
H(w) = w/{l -w) if w < 1/2 and H(w) = (1 - w)/w if 
w > 1/2. But the Lyapunov exponent is still zero. The prob- 
lem now is that the invariant density, p(w) = X/w, is not nor- 
malizable and all the contribution to the Lyapunov exponent 
comes from points that are arbitrarily close to w = 0. Similar 
problems are present in the continuous dynamics since in the 
usual coordinates the model is either nonautonomous or non- 
compact 1 6]. In addition, because the dynamics is limited by a 
cosmological singularity, the domain of frequently used time 
parameters, such as the cosmological time and the volume of 
the universe, is necessarily finite. 

From the above it appears that autonomous equations, 
bounded motions, and normalizable measure are properties 
mutually incompatible in the mixmaster dynamics when all 
the orbits are taken into account. We observe, however, that 
there are invariant bounded subsets of orbits in map F as 
well as invariant subsets with normalizable measure in map 
H which do have positive Lyapunov exponents. For exam- 
ple, map H has a nontrivial set of invariant orbits embedded 
in the interval [a, (1 + a) -1 ], for every a £ (0, 1/3], which 
has normalizable measure. The invariant set is composed of 
all the orbits that never leave this interval and as such con- 
tains a countable number of periodic orbits and an uncount- 
able number of nonperiodic orbits. The Lyapunov exponent 
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of the invariant set, as computed along typical nonperiodic or- 
bits, satisfies h > 21n(l + a). Since map H corresponds to 
the asymptotic behavior of one parametrization of the con- 
tinuous mixmaster model, an invariant set of map H must 
correspond to an invariant set of the asymptotic dynamics for 
this particular parametrization, which therefore satisfies con- 
ditions (z'-z'v) and has positive Lyapunov exponent. It follows 
then from our main result (0 that these asymptotically in- 
variant sets must have positive Lyapunov exponents for any 
space-time reparametrization of the continuous dynamics that 
preserves these conditions on them. These invariant sets are 
therefore chaotic with respect to any coordinate system for 
which the Lyapunov exponents can be properly computed. 
In these sense one can meaningfully say that the mixmaster 
cosmology exhibits coordinate independent chaotic behavior 
close to the big bang. 

Summarizing, we have uncovered the transformation laws 
of the Lyapunov exponents for flows under space-time 
reparametrizations. Strikingly, systems exhibiting exponen- 
tial separation of nearby orbits with respect to one choice of 
the time parameter will display exponential divergence with 
respect to any other time parameter that preserves condi- 
tions (z'-z'v). This implies that chaos is invariant under time 
reparametrizations, which is in sharp contrast with previous 
results in relativistic cosmology, where the apparent nonin- 
variance of chaos has been the subject of an intensive debate. 
Our findings thus shed new light on the conceptual problem 
of chaos in cosmology lE2il . 
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